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The  purpose  of  this  investigation  is  to  analyze  combustion  characteristics  of  biogas  flameless  mode  based 
on  clean  technology  development  strategies.  A  three  dimensional  (3D)  computational  fluid  dynamic 
(CFD)  study  has  been  performed  to  illustrate  various  priorities  of  biogas  flameless  combustion  compared 
to  the  conventional  mode.  The  effects  of  preheated  temperature  and  wall  temperature,  reaction  zone  and 
pollutant  formation  are  observed  and  the  impacts  of  combustion  and  turbulence  models  on  numerical 
results  are  discussed.  Although  preheated  conventional  combustion  could  be  effective  in  terms  of  fuel 
consumption  reduction,  NO*  formation  increases.  It  has  been  found  that  biogas  is  not  eligible  to  be 
applied  in  furnace  heat  up  due  to  its  low  calorific  value  (LCV)  and  it  is  necessary  to  utilize  a  high  calorific 
value  fuel  to  preheat  the  furnace.  The  required  enthalpy  for  biogas  auto-ignition  temperature  is  supplied 
by  enthalpy  of  preheated  oxidizer.  In  biogas  flameless  combustion,  the  mean  temperature  of  the  furnace 
is  lower  than  traditional  combustion  throughout  the  chamber.  Compared  to  the  biogas  flameless  combus¬ 
tion  with  uniform  temperature,  very  high  and  fluctuated  temperatures  are  recorded  in  conventional  com¬ 
bustion.  Since  high  entropy  generation  intensifies  irreversibility,  exergy  loss  is  higher  in  biogas 
conventional  combustion  compared  to  the  biogas  flameless  regime.  Entropy  generation  minimization 
in  flameless  mode  is  attributed  to  the  uniform  temperature  inside  the  chamber. 
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1.  Introduction 

Industrial  development  is  in  debt  of  energy  consumption  and 
more  than  80%  of  the  world  energy  demand  is  supplied  by  different 
kinds  of  fossil  fuels.  However,  the  resources  of  fossil  fuels  are  de¬ 
pleted  day  by  day  and  the  future  energy  scenario  of  the  world 
has  become  one  of  the  main  concerns.  Indeed,  environmental  con¬ 
cerns  have  increased  due  to  raising  rate  of  emissions  released  from 
fossil  fuel  combustion.  Global  warming  (GW)  has  become  one  of 
the  most  important  environmental  issues  due  to  increasing  rate 
of  greenhouse  gases  (GHGs)  generation.  Therefore,  the  request 
for  clean  alternative  fuel  and  efficient  combustion  technology  has 
become  more  important  [1-3].  In  recent  decades,  utilization  of 
renewable  and  sustainable  energy  such  as  biomass,  solar  energy, 
wind  energy,  hydropower  and  geothermal  has  been  developed 
properly.  Furthermore,  biogas  from  wastewater  effluent,  municipal 
solid  wastes  (MSW),  animal  waste  and  agricultural  by-products 
have  been  employed  for  combined  heat  and  power  (CHP)  genera¬ 
tion  purposes  [4-8  .  In  the  other  hand,  necessity  of  biogas  collec¬ 
tion  from  aforementioned  resources  is  unavoidable  because  CH4 
and  C02  as  the  main  components  of  biogas  participate  in  the  GW 
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constitution  actively  [9].  Since  the  negative  effect  of  CH4  on  the 
GW  is  23  times  more  than  C02,  biogas  collection  from  anaerobic 
digestion  (AD)  has  become  more  highlighted.  In  most  of  the  AD, 
the  percentage  of  CH4  is  enough  to  be  considered  as  a  clean  fuel. 
The  amount  of  CH4  in  biogas  components  depends  on  the  feedstock 
is  various  (from  40%  up  to  60%)  [10].  In  biogas  conventional  com¬ 
bustion,  pollutant  formation  mitigates  compared  to  traditional 
combustion  of  pure  CH4.  However,  biogas  traditional  combustion 
encounters  some  problems  due  to  LCV  of  biogas.  Therefore,  biogas 
should  be  upgraded  to  remove  its  non-combustible  C02  impurity 
[11  ].  Based  on  the  application  of  biogas,  pure  biogas  can  be  cleaned 
and  upgraded  by  some  technologies  such  as  water  scrubbing,  cryo¬ 
genic  process  and  membrane.  In  order  to  prevent  implementation 
of  upgrading  equipment,  application  of  flameless  combustion  was 
proposed  for  pure  biogas  combustion.  Since  combustion  is  still  the 
most  important  technique  for  energy  generation,  flameless  com¬ 
bustion  was  introduced  to  improve  the  combustion  efficiency 
and  decrease  pollutant  formation  concomitantly  [12,13  .  The 
importance  of  flameless  combustion  technology  has  become  more 
highlighted  when  the  inability  of  other  combustion  technologies  in 
terms  of  simultaneous  pollutant  reduction  and  thermal  efficiency 
enhancement  was  proven.  Low  emission  formation  in  flameless 
combustion  is  obtained  due  to  dilution  of  the  air  combustion  by 
inert  gases  such  as  nitrogen  (N2)  and  C02  [14-16].  Therefore, 
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Nomenclature 

r 

position  vector 

Ts,o 

outer  surface  temperature 

s 

direction  vector 

Too 

the  ambient  temperature  set  at  300°I< 

I 

radiation  intensity 

h 

natural  convection  coefficient 

h 

radiation  intensity  for  wavelength  2 

G 

Stephane-Boltzmann  constant  (5.67  x  10”8  W/m2  K4) 

spectral  absorption  coefficient 

S 

solid  surface  emissivity 

^bX 

black  body  intensity 

X 

percentage  of  excess  air 

s' 

scattering  direction  vector 

P 

mole  fraction  of  C02  in  biogas 

s 

path  length 

r 

density  (Kg/m3) 

n 

refractive  index 

kf 

thermal  conductivity  of  fluid 

scattering  coefficient 

Yi 

mass  fraction  of  species  i 

G 

Stephan-Boltzmann  constant  (5.669  x  10“8  W/m2  I<4) 

Wi 

generation  or  consumption  of  specious 

0 

phase  function 

Q' 

solid  angle 

combustion  instability  can  be  appeared  due  to  low  amounts  of  oxi¬ 
dizer  and  consequently  flame  quenching  occurs.  Thus  the  high 
temperature  of  diluted  oxidizer  plays  an  important  role  to  exceed 
the  auto  -  ignition  temperature  of  the  fuel.  It  can  be  construed  that 
required  enthalpy  is  supplied  by  preheating  the  oxidizer  to  achieve 
the  self-ignition  temperature  of  the  fuel  [17].  Although  the  great 
performance  of  flameless  combustion  has  been  developed  in  terms 
of  fossil  fuel  utilization,  various  aspects  of  biogas  flameless  com¬ 
bustion  have  not  been  investigated  properly.  Since  computational 
analyses  are  becoming  more  important  due  to  their  acceptable 
accuracy  and  lower  cost,  different  perspectives  of  biogas  flameless 
combustion  can  be  analyzed  numerically.  The  theoretical  study 
using  CFD  software  decreases  errors  and  trials  on  experimental 
investigations  for  development  of  new  models.  In  this  article,  a 
CFD  modeling  of  biogas  flameless  combustion  is  investigated  to 
identify  the  most  important  parameters  that  should  be  considered 
in  biogas  flameless  regime.  The  effects  of  chemical  kinetic  mecha¬ 
nisms,  turbulence  models,  radiation  heat  transfer  and  differential 
diffusion  effects  on  the  accuracy  of  the  model  are  noted. 

2.  Conventional  combustion  of  LCV  fuels 

One  of  the  most  important  barriers  for  biogas  utilization  devel¬ 
opment  is  its  low  calorific  value  as  well  as  its  corrosive  nature.  LCV 
fuels  were  not  taken  into  account  in  the  energy  mix  of  the  world 
due  to  the  abundance  of  high  calorific  value  fossil  fuels  and  lower 
energy  prices  in  the  past.  However,  recent  enhancements  of  energy 
prices  and  concerns  of  the  future  energy  scenario  have  attracted 
more  attention  to  LCV  fuel  utilization  in  power  generation.  There¬ 
fore,  combustion  characteristics  of  various  combinations  of  biogas 
have  become  more  important.  It  has  been  claimed  that  raw  biogas 
cannot  be  utilized  in  commercial  burners  directly  as  the  substitu¬ 
tion  of  natural  gas  (NG)  or  liquid  petroleum  gas  (LPG).  If  biogas 
conducted  to  the  traditional  burner  orifice  at  the  pressure  intended 
for  charging  NG  or  LPG,  the  air  fuel  ratio  is  not  sufficient  for  com¬ 
bustion  stability  due  to  high  C02  content  of  biogas.  Thus,  a  new 
burner  with  the  separate  control  measurement  system  should  be 
installed  for  biogas.  Since  water  vapor  and  H2S  as  the  components 
of  biogas  have  corrosive  characteristics,  the  condensation  should 
be  prevented  by  maintaining  the  furnace  temperatures  above  the 
dew  point  temperature.  It  means  that  the  combustion  furnace 
should  be  heated  up  by  NG  or  LPG  to  obtain  higher  operating  tem¬ 
perature  [18].  Therefore,  a  dual  role  burner  should  be  installed  in 
conventional  combustion  system  fueled  by  biogas.  Indeed,  in  order 
to  switch  over  from  NG  or  LPG  to  biogas  some  equipment  should 
be  employed.  Consequently,  complicated  setting  and  the  low 
efficiency  of  the  conventional  combustion  systems  can  encourage 


biogas  users  to  apply  flameless  combustion  in  their  combustion 
process  [19]. 

3.  Numerical  studies  of  flameless  regime 

Flameless  combustion  has  been  referred  as  flameless  oxidation, 
colorless  distributed  combustion,  moderate  or  intense  low  oxygen 
dilution  (MILD)  combustion  and  high-temperature  air  combustion 
(HiTac)  [20].  In  the  flameless  mode,  combustion  takes  place  in  a 
distributed  reaction  zone  with  uniform  low  temperature.  The  fluc¬ 
tuations  of  temperature  are  omitted  in  flameless  regime  compared 
to  the  conventional  flames.  Combustion  occurs  in  a  low  oxygen 
concentration  atmosphere  without  visible  flame.  Indeed,  the  levels 
of  the  noise,  NOx  and  soot  emissions  decrease  [21  .  Due  to  compli¬ 
cated  conditions  in  flameless  mode,  numerical  modeling  of  the 
flameless  regime  has  received  especial  attention.  Weber  et  al. 
[22]  modeled  an  industrial  chamber  with  square  cross  section  in 
steady  state  flameless  combustion  condition.  Very  high  momen¬ 
tum  oxidizer  was  charged  through  a  central  hole,  NG  was  entered 
to  the  furnace  through  two  injectors  as  a  fuel  and  exhaust  gases 
were  recirculated  in  the  system.  In  these  circumstances,  turbu¬ 
lence  was  modeled  by  standard  k-s  model  [23,24]  and  RNG  k-s 
model  [25].  The  combustion  model  was  simulated  with  the  eddy 
breakup  model  and  a  two-step  reaction  scheme;  the  eddy  dissipa¬ 
tion  model  (EDM)  with  the  chemical  equilibrium  and  the  con¬ 
served  scalar/prescribed  probability  density  function  (pdf)  with  a 
chemical  equilibrium  assumption.  A  similar  numerical  investiga¬ 
tion  was  done  by  Kim  et  al.  [26]  who  applied  the  standard  k-s 
model,  the  eddy-dissipation  concept  (EDC)  and  four  different  glo¬ 
bal  reaction  mechanisms.  The  characteristics  of  the  MILD  combus¬ 
tion  system  included  a  burner  with  a  central  fuel  jet  surrounded  by 
six  oxidizer  jets  and  exhaust  gases  recirculation  system  was  inves¬ 
tigated  numerically  using  k-s  and  flamelet  model  by  Coelho  and 
Peters  et  al.  [27]  and  Dally  et  al.  [28].  It  was  pointed  out  that  the 
numerical  temperatures  along  the  centerline  of  combustor  were 
in  good  agreement  with  experimental  records  for  pure  CH4  and 
CH4  diluted  with  C02,  and  in  poor  agreement  for  CH4  diluted  with 
N2.  Khoshhal  et  al.  [29]  investigated  NOx  formation  and  heat  trans¬ 
fer  mechanisms  in  a  high  temperature  air  combustion  boiler 
numerically.  It  was  pointed  out  that  the  experimental  measured 
values  and  the  CFD  results  illustrated  good  agreement.  Although 
the  temperatures  of  preheated  oxidizer  in  flameless  combustion 
mode  are  higher  than  conventional  conditions,  the  reaction-con- 
trolling  temperatures  are  lower  than  traditional  combustion  due 
to  the  low  oxygen  concentration.  In  the  other  hand,  burner  config¬ 
uration  is  very  important  to  achieve  flameless  conditions.  In  some 
investigations,  burner  contains  a  central  inlet  jet  of  fuel  and  some 
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inlet  preheated  air  jets  located  around  the  central  jet  [30].  Also,  the 
high-momentum  preheated  oxidizer  can  be  conducted  from  the 
central  of  the  burner  surrounded  by  some  low-momentum  fuel  jets 
[31  .  The  characteristics  of  the  turbulent  non-premixed  methane/ 
hydrogen  flame  issued  by  a  jet  in  hot  co-flow  (JHC)  were  investi¬ 
gated  numerically  by  Christo  and  Dally  [32  .  It  was  stipulated  that 
differential  diffusion  effects  should  be  considered  in  CFD  modeling 
due  to  their  crucial  role  in  the  accuracy  of  the  calculations.  Indeed, 
the  numerical  results  of  detailed  chemical  kinetics  like  the  eddy 
dissipation  concept  (EDC)  solver,  illustrated  higher  accuracy  in 
JHC.  Galleti  et  al.  [33]  simulated  a  recuperative  burner  where  the 
fuel  was  charged  at  the  burner  axis  and  the  oxidizer  was  injected 
in  a  co-flow  configuration.  A  flame  tube  with  three  windows  was 
set  to  promote  the  internal  recirculation  of  exhaust  gases.  Exhaust 
gas  recirculation  was  done  through  those  windows  which  mixed 
with  the  air  before  combustion.  A  modified  standard  k-s  model, 
a  single  one  step  mechanism  and  finite  rate  eddy  dissipation  meth¬ 
od  were  considered  in  this  study.  Indeed,  Parente  et  al.  [34]  studied 
NO*  formation  in  the  aforementioned  combustion  system  numeri¬ 
cally.  Various  aspects  of  turbulence  flameless  models  in  high  pres¬ 
sure  gas  turbine  were  investigated  by  Schiitz  et  al.  [35  .  The 
premixed  oxidizer  and  fuel  were  fed  into  the  chamber  via  12  noz¬ 
zles  located  along  a  circle.  Thanks  to  the  mass  continuity  and  very 
high  velocity  of  the  mixture,  a  recirculating  annulus  was  formed 
around  the  nozzles.  CFD  simulation  of  the  system  was  carried 
out  either  by  GRI3.0  or  by  two-step  global  mechanism  using  k-s 
method  as  the  viscous  model  and  EDC  as  the  turbulence  chemistry 
interaction  model.  Although,  the  numerical  results  were  in  good 
agreement  with  experimental  records  some  shortcomings  in  the 
estimation  of  turbulence  mixing  were  reported.  In  the  other  at¬ 
tempt,  the  characteristics  of  LCV  gases  combustion  in  a  similar  fur¬ 
nace  operated  at  atmospheric  temperature  was  simulated  by 
Danon  et  al.  [36].  Both  realizable  k-s  model  and  the  Reynolds 
stress  turbulence  were  applied  in  this  modeling.  An  EDC  model 
with  two  various  chemical  mechanisms  were  set.  The  differences 
between  mean  temperatures  along  the  axis  in  numerical  and 
experimental  results  were  reported  in  the  uncertainty  range  of 
the  optimal.  Chen  and  Zheng  [37]  investigated  the  effects  of  oxy¬ 
gen  concentration,  the  temperature  of  oxidizer  and  the  percentage 
of  hydrogen  on  the  hydrogen-enriched  biogas  counter-flow  flame¬ 
less  combustion  mathematically.  It  was  pointed  out  hydrogen-en¬ 
riched  biogas  flameless  combustion  can  be  sustained  with  low 
preheated  temperature  of  the  oxidizer  and  extremely  highly  di¬ 
luted  oxygen  concentration  in  the  oxidizer  mixture.  They  applied 
flamelet  approach  for  modeling  the  combustion  and  found  an  ur¬ 
gency  to  develop  new  turbulent  combustion  methods  for  flameless 
combustion. 

Based  on  aforementioned  literatures,  it  can  be  construed  that 
the  numerical  simulation  of  the  flameless  combustion  method  is 
a  challenging  problem,  due  to  the  lack  of  comprehensive  theoreti¬ 
cal  knowledge  about  this  technique.  Although  excellent  agree¬ 
ments  between  experimental  and  numerical  investigations  have 
sometimes  been  reported,  few  investigations  were  presented  a 
holistic  comparison  of  various  turbulence  and  combustion  models, 
chemical  reaction  mechanisms,  and  compared  their  numerical  re¬ 
sults  with  experimental  data.  Furthermore,  a  lack  of  comprehen¬ 
sive  numerical  investigation  about  biogas  flameless  combustion 
is  observed  due  to  complicated  chemical  reaction  mechanism. 
Though  the  concept  of  fossil  fuel  flameless  combustion  has  been 
extensively  investigated  numerically  and  mathematically,  biogas 
flameless  modeling  has  received  little  attention.  Since  low  oxygen 
concentration  of  biogas  flameless  combustion  has  led  the  system 
to  slower  reaction  rates  and  increase  the  effects  of  molecular  diffu¬ 
sion  on  combustion  characteristics,  the  applications  of  combustion 
models  that  presume  fast  chemistry  and  eliminate  the  impacts  of 
differential  diffusion  are  challenged  [38  . 


4.  CFD  modeling 

The  CFD  model  of  combustor  designed  for  this  numerical  study 
is  based  on  the  geometry  of  the  experimental  lab-scale  biogas 
flameless  combustion  carried  out  by  Hosseini  et  al.  [39  .  The  inside 
diameter  and  the  length  of  the  chamber  are  150  mm  and  600  mm 
respectively.  The  5  mm  central  inlet  of  the  burner  is  considered  for 
biogas  inlet  and  the  other  inlets  are  the  oxidizer  entrances.  The  ex¬ 
haust  gases  are  conducted  to  the  outside  of  the  chamber  through  a 
central  hole  with  50  mm  diameter.  Fig.  1  depicts  the  schematic  of 
the  biogas  flameless  furnace. 

This  3D  simulation  is  performed  with  ANSYS  14  using  ANSYS 
Modeler  to  design  the  flameless  furnace  and  ANSYS  Meshing  to 
mesh  the  furnace  [40].  Convergence  rate  as  well  as  scalar  proper¬ 
ties  can  be  improved  by  mesh  refinement  and  grid  resolution  for 
smooth  flow  representation  can  be  ensured.  The  geometry  of  the 
furnace  is  not  complicated  in  this  numerical  furnace  model,  there¬ 
fore  calculations  speed  increase  due  to  mitigate  meshing  nodes 
and  elements.  The  quantity  of  mesh  grids  directly  impacts  on  the 
solution  duration.  Since  the  chamber  model  is  symmetrical,  just 
an  eighth  part  of  the  model  is  solved.  Fig.  2  demonstrates  the 
meshes  of  the  flameless  chamber.  To  enhance  the  precision  of  pre¬ 
dictions,  the  regions  near  the  air  and  fuel  nozzle  have  smaller  con¬ 
trol  volume  meshes. 

The  results  of  simulation  with  7769  nodes  and  33,798  elements 
is  in  good  agreement  with  the  experimental  records.  The  grid  inde¬ 
pendent  of  the  simulation  was  tested  by  changing  the  number  of 
nodes  to  the  finer  meshes.  The  number  of  elements  were  adopted 
33,798,  68,453  and  100,356  to  find  the  most  conformity  to  the 
experimental  records.  However  meaningful  changes  were  not 
observed  in  the  results  and  the  grid  independent  of  solution  is 


Fig.  1.  Schematic  of  the  flameless  furnace. 


Fig.  2.  An  eighth  part  of  the  flameless  furnace  after  meshing. 
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Axial  distance  [m] 

Fig.  3.  The  axial  temperature  variation  of  experimental  and  numerical  results. 


confirmed  by  this  test.  Fig.  3  shows  the  axial  temperature 
variations  of  experimental  and  simulated  data  with  respect  to 
the  various  meshes  when  preheated  air  temperature  was  900  K 
and  the  concentration  of  oxygen  in  the  oxidizer  was  set  7%.  This 
figure  confirms  that  the  experimental  and  simulated  data  have 
the  same  trend  in  terms  of  the  temperature  profile,  and  mesh  with 
7769  nodes  was  selected  due  to  lower  computational  cost. 

The  governing  equations  (transport  equations  include  continu¬ 
ity,  energy,  momentum)  are  solved  using  CFD  package  ANSYS  FLU- 
ENT14  [40].  In  Reynolds  averaging,  the  variables  in  the  exact 
Navier-Stoks  equations  are  converted  into  he  mean  and  fluctuating 
components.  For  all  scalar  quantities  like  energy,  pressure  and 
species  concentration  it  can  be  written: 

0  =  0  +  0'  (1) 

Likewise,  for  velocity: 

Ui  =  Ui  +  u[  (2) 


where  u,  and  u j  are  mean  and  fluctuating  velocity  respectively  (i  =  1, 
2,3). 

The  continuity  and  momentum  equations  called  Reynolds-aver¬ 
age  Navier-Stokes  (RANS)  in  turbulent  circumstance  are  written  as 
follows  equations  respectively. 


dp 

~dt 


+ 


d_ 

dXi 


(pUi)  =  o 


d  d  dp  d 


dt 


+ 


dXi  dXj 

d 

dXj 


dllj  dUj  2  d  llf 

dXi  +  dXj  3  ij  dXi 


(-pu'iU'j) 


(4) 


Also,  energy  transfer  due  to  conduction,  species  diffusion  and 
viscous  dissipation  is: 


|  ( pE )  +  V  •  (v(pE  +  p))  =  V  •  L keff  vr  -  Y)hJj  +  ( xeff  •  V)  j  +  s„ 

(5) 

where  /+//  is  effective  conductivity  and])  is  diffusion  flux  of  species  j. 

On  the  right  hand  side  of  Eq.  (5),  the  first  three  terms  represent 
energy  transfer  due  to  conduction,  species  diffusion  and  viscous 
dissipation  respectively  and  Sh  is  a  combination  of  heat  chemical 
reactant  and  every  other  heat  resource  which  are  defined. 

Second-order  discretization  scheme  is  applied  to  solve  all  gov¬ 
erning  equations.  The  residual  of  the  energy  equation  should  drop 


below  10-6  and  for  all  other  variables  it  is  set  at  10-3  to  ensure  the 
convergence  of  the  solution.  The  swirl  velocity  of  components  is 
eliminated  in  this  steady-state  CFD  study.  To  investigate  various 
aspects  of  biogas  flameless  combustion,  the  simulation  is  done  in 
stoichiometric  equivalence  ratios  =  1)  and  the  results  are  com¬ 
pared  by  conventional  combustion.  Equivalence  ratio  is  applied 
to  show  quantitatively  whether  the  fuel-air  in  a  chemical  reaction 
is  lean  ( <P  <  1 ),  stoichiometric  ( <3?  =  1 )  or  rich  ( <P  >  1 )  [41  .  The  den¬ 
sity  of  biogas  consists  of  60%  CH4  and  40%  C02  in  300  K  is  consid¬ 
ered  1.106  Kg/m3.  Table  1  demonstrates  the  densities  of  preheated 
air  for  different  temperatures  and  various  mixtures  of  N2  and  02 
The  global  combustion  reaction  of  biogas,  based  on  a  different 
mole  fraction  of  C02  is  presented  in  Eq.  (6)  [42  . 

(1  -  P)[CH4  +  (2  +  x)(02  +  3.76N2)]  +  pC02 
— >  C02  +  2(1  -  P)  +  H20  +  7.56(1  -  /?)N2  +  x(l  -  P)(02 

+  3.76N2)  (6) 

Which  x  is  the  percentage  of  excess  air  and  p  is  the  mole  frac¬ 
tion  of  C02  in  biogas. 

For  biogas  consist  of  40%  C02  and  60%  CFI4,  Eq.  (2)  can  be  writ¬ 
ten  as  Eq.  (7). 

0.6CH4  +  0.4CO2  +  1 .2(02  +  3.76N2) 

— >  C02  +  1.2H20  +  4.512N2  (7) 

Based  on  a  two  step  model  developed  by  Westbrook  and  Dryer 
et  al.  [43  (WD  model)  for  CH4  combustion,  the  model  consists  of 
two  following  chemical  reactions. 

CH4  +  3/202  -►  CO  +  2H20  (8) 

CO  +  l/202  C02  (9) 

In  the  chemical  reaction  (8),  carbon  monoxide  (CO)  and  water 
vapor  (H20)  are  formed,  while  in  the  chemical  reaction  (9)  oxida¬ 
tion  of  CO  to  C02  and  its  dissociation  takes  place.  This  WD  reduced 
model  illustrates  CO,  C02  and  F120  formation.  A  chemical  mecha¬ 
nism  for  combustion  chemistry  modeling  has  at  least  three  parts: 
gas  phase  kinetic  file  (a  list  of  all  chemical  reactions  of  interest 
to  simulation  including  appropriate  Arrhenius),  thermodynamic 
data  base  (consist  of  thermodynamic  coefficients  for  all  gas  phase 
kinetic  file)  and  transport  data  file.  The  coefficient  used  for  Eqs.  (8) 
and  (9)  are  values  of  Arrhenius  equation  coefficients. 

kfi  =  A,TBi  exp  (10) 

where  A,  is  pre-exponential  factor  (for  Eq.  (8),  it  is  equal  to  2.8E+9 
and  for  Eq.  (9)  it  is  10.00E+14),  pi  is  temperature  coefficient  (equal 
zero  for  both  Eqs.  (8)  and  (9))  and  E,  is  activation  energy  (for  Eq.  (8), 
it  is  equal  to  48.4  and  for  Eq.  (9)  it  is  40). 


Table  1 

The  density  (kg/m3)  of  diluted  air  in  various  temperatures. 


Temperature  (K) 

5%02 

7%02 

10%O2 

15%02 

21%02 

300 

1.146 

1.15 

1.154 

1.162 

1.177 

400 

0.859 

0.862 

0.866 

0.8717 

0.8829 

500 

0.687 

0.689 

0.692 

0.6973 

0.7063 

600 

0.573 

0.575 

0.577 

0.5811 

0.5886 

700 

0.491 

0.492 

0.4946 

0.4981 

0.5046 

800 

0.429 

0.431 

0.433 

0.4358 

0.4415 

900 

0.382 

0.383 

0.3847 

0.3874 

0.3925 

1000 

0.344 

0.345 

0.3462 

0.3486 

0.3532 

1100 

0.312 

0.313 

0.3147 

0.3169 

0.3211 

1200 

0.286 

0.287 

0.2885 

0.2905 

0.2944 

S.E.  Hosseini  et  al.f  Energy  Conversion  and  Management  81  (2014)  41-50 


45 


Presented  gas  phase  kinetic  file  with  thermodynamic  data  base 
and  transport  data  file  were  applied  for  further  simulation  as  re¬ 
duced  WD  chemical  reaction  mechanism. 

The  Reynolds-averaged  from  of  the  governing  equation  is  ap¬ 
plied  in  conjunction  with  a  turbulence  model.  For  turbulent  kinetic 
energy  (/<)  and  its  dissipation  (s),  the  standard  k-s  formulation  is 
considered  because  of  its  accuracy  and  robustness  for  a  wide  range 
of  turbulent  flows.  Indeed,  since  the  k-s  model  is  valid  for  fully  tur¬ 
bulent  flow  in  CFD  modeling  and  turbulence  condition  is  effective 
in  flameless  combustion  formation  inside  the  chamber,  this  model 
could  increase  the  accuracy  of  the  calculations  [44].  The  thermal 
conductivity,  specific  heat  and  viscosity  are  computed  with  respect 
to  the  average  of  species  mass  fraction.  Since  the  length  of  the 
chamber  is  longer  than  the  molecular  mean-free  path  of  the  reac¬ 
tants  flowing  through  the  chamber,  Navier-Stokes  equations  are 
valid  in  the  simulation.  Radiation  heat  transfer  between  furnace  in¬ 
ner  walls  is  taken  into  account.  A  radiative  transfer  equation  is 
solved  for  discrete  solid  angles  across  the  numerical  domain.  Also, 
spatial  variation  in  the  total  emissivity  is  calculated  as  a  function  of 
biogas  composition  and  temperature  because  the  weighted  sum  of 
the  gray  gas  model  (WSGGM)  is  incorporated.  In  WSGGM  the  total 
emissivity  is  considered  as  a  function  of  the  H20  and  C02  temper¬ 
ature  and  local  mass  fractions.  The  heat  transfer  radiation  for  spec¬ 
tral  intensity  for  position  r  in  direction  s(/;  (r,s))  is: 

V  •  (h(r,s),s)  +  (ax  +  as)h(r,s) 
rr  r4n 

=  a,n2Ibi  +  -^_  /  r(f,s')<b(s,s')dQ!  (11) 

0 

Since  simulation  of  radiation  heat  transfer  from  various  thick¬ 
nesses  is  possible  by  discrete  ordinates  (DO)  model,  surface-to-sur- 
face  radiation  heat  transfer  in  the  studied  flameless  combustor  is 
modeled  by  DO.  At  the  internal  walls,  non-slip  and  no  species  flux 
normal  to  the  surfaces  are  considered.  Moreover,  heat  loss  from 
combustor  surfaces  to  the  surroundings  is  taken  into  account 
through  Eq.  (12)  in  which,  both  thermal  heat  transfer  radiation 
and  natural  convection  are  accounted. 

<7  =  h(Ts,o  ~  Too)  +  £(T(T4so  -  O  (12) 

where  Ts  o  is  the  outer  surface  temperature,  T ^  is  the  ambient  tem¬ 
perature  set  at  300  K,  h  is  natural  convection  coefficient  considered 
constant  value  5  W/m2  K,  a  =  5.67  x  1  O'8  W/m2  K4  is  Stephane- 
Boltzmann  constant  and  s  is  the  solid  surface  emissivity. 

The  velocity  distribution  of  the  inlet  oxidizer  and  fuel  are  uni¬ 
form.  The  temperature  of  the  inlet  biogas  is  300  K  and  the  pressure 
outlet  is  1.013  x  105  Pa.  Eddy  dissipation  model  (EDM)  and  volu¬ 
metric  reaction-based  with  respect  to  the  two-step  chemical  ki¬ 
netic  mechanisms  of  biogas  combustion  is  set.  Therefore  the 
application  of  an  ignition  source  is  not  necessary  to  ignite  the  reac¬ 
tants  and  combustion  proceeds  in  the  presence  of  turbulence  con¬ 
dition  (k/s  >  0).  In  the  literature  review  it  was  stipulated  that 
simulation  of  flameless  combustion  with  EDC  shows  better  results 
in  the  previous  documents.  Flowever,  in  some  cases,  it  has  been 
mentioned  that  EDC  is  not  appropriate  for  simulation.  For  instance 
Cristo  and  Dally  et  al.  [32  pointed  out  that  in  some  simulated 
cases  EDC  has  overestimation.  Damkohler  number  distribution 
(Da),  which  represents  the  flow  to  chemical  time-scale  ratio  could 
be  effective  for  combustion  model  selection.  Large  Damkohler  va¬ 
lue  indicates  mixing  controlled  flames  and  low  Damkohler  value 
corresponds  to  slow  chemical  reactions.  It  means  that  reactants 
and  products  are  quickly  combined  by  turbulence  and  the  combus¬ 
tion  system  behaves  like  a  perfect  stirred  reactor.  Therefore,  calcu¬ 
lation  of  the  Damkohler  number  is  required  to  set  proper  flow  and 
chemical  time-scales.  The  Damkohler  (Da)  number  characterizes 
the  behavior  between  mixing  and  reaction  in  a  system,  given  by 
the  ratio  of  a  mixing  or  flow  time-scale  to  a  chemical  time-scale 


(t//tc)  [45  .  Based  Ref.  [40],  EDM  is  appropriate  when  Da  >  1.  Since 
in  this  study  the  distance  between  the  centerline  of  fuel  and  oxi¬ 
dizer  inlet  is  52.5  mm,  mixing  time  is  long  and  chemical  reaction 
time  is  short  (due  to  biogas  utilization)  and  Da  is  greater  than 
one.  Thus,  application  of  EDM  shows  acceptable  results  in  this  spe¬ 
cific  geometry.  The  radial  inlet  velocity  of  the  oxidizer  and  fuel  is 
negligible  at  the  inlet  of  the  furnace  boundaries  and  the  uniform 
axial  velocities  of  biogas  and  oxidizer  are  set  from  the  calculated 
biogas  and  oxidizer  mass  with  respect  to  the  density  of  the  compo¬ 
nents.  Turbulence  intensity  which  is  defined  as  Eq.  (13)  is  set  10% 
and  20%  for  oxidizer  and  fuel  respectively. 


where  u!  is  the  root-mean-square  of  the  turbulent  velocity  fluctuations 
and  U  is  the  mean  velocity  in  averaged  Reynolds.  In  high  turbulence 
models,  turbulence  intensity  is  considered  5-20%  46  .  The  dissipation 
rate  inlet  values  of  the  turbulent  kinetic  energy  are  set  based  on  the 
Versteeg  and  Malalasekera  47  applying  the  inner  diameter  of  the  bio¬ 
gas  inlet  and  the  hydraulic  diameter  of  the  oxidizer  nozzle  as  charac¬ 
teristic  lengths.  For  pressure-velocity  coupling,  coupled  scheme  is 
applied  in  solution  methods  and  for  other  special  discretization,  sec¬ 
ond  order  upwind  scheme  is  set.  NO*  formation  mechanisms  include 
thermal  NO*  and  prompt  NO*  were  adopted  based  on  the  Ansys  Fluent 
defaults.  The  transport  equation  for  NO  species  formation  is: 

^  ( pYm )  +  V  •  (pvYm)  =  V  •  (pDVYN0)  +  SN0  (14) 

ot 

SN0  includes  thermal  NOx  determined  by  Zeldovich  equations 
and  prompt  NOx.  Zeldovich  equations  are  [48]: 

0  +  N2^NO  +  N  (15) 

N  +  02  NO  +  O 
N  +  OH  <-►  NO  +  H 

The  reaction  rate  constants  of  aforementioned  reactions  are 
selected  from  Ref.  [49]  and  partial  equilibrium  approach  is  taken 
into  consideration  to  compute  the  concentration  of  OH  and  O 
radicals.  The  calculation  of  prompt  NO*  formation  is  done  from 
global  model  presented  in  Ref.  [50]. 

CH  +  N2^HCN  +  N  (16) 

N  +  02  <->  NO  +  O 

HCN  +  OH  <r~>  CN  +  H20 
CN  +  02  ^  NO  +  CO 

The  computed  time  averaged  NOx  results  is  applied  in  Eq.  (14). 
The  computations  were  done  by  a  core  i7  computer  with  16  GB  of 
RAM  and  the  cases  solved  between  the  parallel  cores. 

5.  Results  and  discussion 

5.2.  Heat  up 

In  order  to  achieve  flameless  mode,  the  furnace  should  be  heated 
up  at  the  first  step  and  inside  temperature  of  the  chamber  should  in¬ 
crease  over  the  auto-ignition  of  biogas.  Fig.  4  displays  the  volumetric 
temperature  distribution  in  conventional  combustion  of  the  pure 
CH4  and  biogas.  Also,  the  variation  of  temperature  on  the  centerline 
of  the  furnace  fueled  by  CH4  and  biogas  is  demonstrated  in  this  fig¬ 
ure.  The  distribution  of  the  temperature  inside  the  chamber  con¬ 
firms  that  biogas  is  not  eligible  to  be  applied  for  heating  up 
process  practically.  In  the  other  word,  biogas  conventional 


46 


S.E.  Hosseini  et  al./ Energy  Conversion  and  Management  81  (2014)  41-50 


Temperature 
Volume  Rendering  1 
2.177e+003 


1  708e+003 


1.238e+003 


7.692e+002 

[K] 


Biogas  conventional  combustion 

(A) 


(B) 


Fig.  4.  (A)  Temperature  distribution  in  conventional  combustion  of  the  CH4  and 
biogas.  (B)  Variation  of  temperature  on  the  centerline  of  the  furnace. 


CH4  conventional  combustion 


combustion  cannot  generate  high  temperatures  inside  the  furnace. 
Therefore,  in  order  to  achieve  flameless  mode,  furnace  should  be 
heated  up  by  a  high  calorific  value  fuel  such  as  CH4. 

5.2.  Biogas  conventional  combustion  with  preheated  air 

The  mean  temperature  of  conventional  preheated  combustion 
increases  inside  the  chamber  due  to  high  enthalpy  of  preheated 
air,  consequently  fuel  consumption  reduces  drastically.  Fig.  5  illus¬ 
trates  the  effects  of  preheated  air  temperature  on  the  fuel  con¬ 
sumption  in  conventional  combustion. 

Although,  preheated  conventional  combustion  could  be  useful 
in  fuel  consumption  minimization,  NO*  formation  increases  in 
highly  preheated  combustion.  This  phenomenon  is  attributed  to 
augmentation  of  the  peak  temperature  in  the  furnace,  because 
thermal  NO*  intensifies  sharply  in  high  temperatures  [17  .  Fig.  6 
shows  the  variation  of  peak  temperature  inside  the  furnace  and 
NO*  formation  in  conventional  preheated  combustion  with  respect 
to  the  preheated  oxidizer  temperature.  From  this  figure,  it  can  be 
concluded  that  preheated  conventional  combustion  of  biogas  can¬ 
not  increase  the  inside  temperature  of  the  furnace  significantly. 
However,  conspicuous  variation  can  be  seen  in  the  maximum  tem¬ 
perature  and  NOx  formation  in  CH4  preheated  conventional 
combustion. 
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Fig.  5.  The  effects  of  preheated  air  temperature  on  the  fuel  consumption  in 
conventional  combustion. 
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Fig.  6.  (A)  The  variation  of  peak  temperature  in  conventional  preheated  combus¬ 
tion.  (B)  NOx  formation  in  conventional  preheated  combustion. 


5.3.  Flameless  combustion 

In  this  simulation,  the  concentration  of  oxygen  in  oxidizer  is 
considered  5%,  7%,  10%  and  15%  by  volume  respectively  and  it  is 


assumed  that  oxidizer  is  diluted  by  N2.  Fig.  7  demonstrates 
the  flameless  constitution  in  the  furnace  when  Viniet-air  =  30  m/s, 
^inlet-oxidizer  —  900  K  and  0  —  1 . 

Our  experimental  observations  confirm  that  a  small  flame  is 
constitute  when  oxygen  concentration  in  oxidizer  increased  upto 
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Fig.  7.  Biogas  flameless  combustion  formation  with  respect  to  the  various  oxygen 
concentrations. 


15%.  The  maximum  enthalpy  in  biogas  flameless  combustion  is  re¬ 
corded  648,253  J/Kg  but  in  conventional  combustion  the  maximum 
enthalpy  is  48,089  J/Kg.  It  can  be  construed  that  required  enthalpy 
for  biogas  auto-ignition  temperature  is  supplied  by  preheated  oxi¬ 
dizer.  In  biogas  flameless  combustion,  the  mean  temperature  of  the 
furnace  is  lower  than  CH4  traditional  combustion  throughout  the 
chamber.  Compared  to  the  biogas  flameless  combustion  with  uni¬ 
form  temperature,  very  high  and  fluctuated  temperatures  are  re¬ 
corded  in  conventional  combustion.  Consequently,  hot  spots 
which  play  crucial  role  in  NO*  formation  are  eliminated  in  flame¬ 
less  mode.  Since  C02  constitutes  40%  of  biogas  components,  the 
temperature  of  the  reactants  decreases.  It  can  be  attributed  to 
C02  high  heat  capacity  especially  in  high  temperatures  [39  .  The 
temperature  inside  the  biogas  flameless  furnace  is  uniform  aver¬ 
agely  1090  K.  Also,  wall  temperature  is  lower  and  more  uniform 
in  flameless  regime  compared  to  the  conventional  mode.  Fig.  8 
shows  the  distribution  of  wall  temperature  in  conventional  and 
flameless  regimes. 

Based  on  mathematic  formula  of  temperature  uniformity  ratio 
defined  by  Yang  et  al.  [51],  the  ratio  of  uniformity  in  the  flameless 
furnace  is  calculated  around  zero.  This  condition  not  only  increases 
the  durability  of  the  refractory  inside  the  chamber  but  also  pro¬ 
vides  an  ideal  circumstance  for  NOx  formation  reduction.  Further¬ 
more,  fuel  consumption  decreases  from  3.24  g/s  in  biogas 
conventional  combustion  to  1.07  g/s  in  flameless  mode.  Fig.  9  com¬ 
pared  the  entropy  generation  in  centerline  of  the  furnace  in  tradi¬ 
tional  and  flameless  combustion.  Since  high  entropy  generation 
intensifies  irreversibility,  exergy  loss  is  higher  in  conventional 
combustion.  Som  and  Datta  [52]  pointed  out  that  the  most  irre¬ 
versibility  in  a  combustion  system  is  related  to  the  internal  heat 
transfer  within  the  combustor  between  the  products  and  reactants. 
Chen  et  al.  [53]  stipulated  that  in  the  high  temperature  air  combus¬ 
tion  the  irreversibility  is  attributed  to  the  heat  transfer  in  the  reac¬ 
tion  zones,  however  in  the  flameless  mode,  the  reaction  zones  are 
dominated  by  the  irreversibility  due  to  chemical  reaction.  Entropy 
generation  minimization  in  flameless  mode  could  be  attributed  to 
the  uniform  temperature  inside  the  chamber. 
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Fig.  8.  Wall  temperature  distribution. 
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Fig.  9.  Entropy  generation  in  the  furnace. 


5.4.  Species  concentration  in  biogas  flameless  chamber 

Fig.  10  demonstrates  concentration  of  CH4  and  C02  along  the 
central  axis  of  the  biogas  flameless  furnace.  Since  CH4  and  C02 
form  60%  and  40%  of  biogas  respectively,  the  concentration  of  these 
species  is  very  high  at  the  entrance  of  the  chamber.  However,  the 
concentrations  of  CH4  and  C02  decrease  sharply  in  biogas  flameless 
mode  compared  to  the  conventional  combustion.  It  means  that 
combustion  phenomena  takes  place  thorough  the  furnace  homoge¬ 
neously  in  flameless  mode. 

The  streamline  of  the  species  in  both  conventional  and  flame¬ 
less  mode  is  shown  in  Fig.  11.  The  species  of  biogas  are  pulled  to 
the  oxidizer  streamline  due  to  various  pressures.  Since  the  velocity 
of  the  oxidizer  is  more  than  biogas,  the  vacuum  circumstance 
causes  the  biogas  species  to  be  mixed  with  oxidizer.  The  mixture 
of  the  fuel  and  oxidizer  species  in  flameless  regime  is  done  sooner 
than  conventional  mode  due  to  the  turbulence  environment.  In  the 
other  word,  the  complete  biogas  combustion  occurred  in  the  zone 
close  to  the  burner. 

Fig.  12  confirms  that  the  concentration  of  02  along  the  furnace 
is  very  low  and  uniform.  The  low  concentration  of  02  can  be  attrib¬ 
uted  to  diluted  oxidizer  and  the  complete  and  quick  mixing  pro¬ 
cess  in  biogas  flameless  mode.  Apart  from  that,  the  uniform 
temperature  inside  the  flameless  chamber  which  is  higher  than 
the  biogas  self-ignition  temperature  guarantees  the  stable  condi¬ 
tions  for  the  biogas  flameless  regime  formation.  The  high  amount 
of  CO  concentration  in  flameless  regime  is  attributed  to  low  OH 
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Fig.  10.  Concentration  of  CH4  and  C02  along  the  central  axis  of  the  furnace. 
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Fig.  12.  (A)  Oxygen  concentration  and  (B)  monoxide  carbon  concentration  pattern 
inside  the  furnace. 


radical  concentration  which  controls  the  CO  conversion  to  C02 

[54], 

Fig.  13  depicts  NO*  formation  pattern  along  the  chamber.  The 
constitution  of  NOx  is  very  low  in  flameless  mode.  However,  in  con¬ 
ventional  combustion  of  CH4  and  biogas  the  NO*  emission  in  the 
outlet  of  the  chamber  was  recorded  45  ppm  and  22  ppm  respec¬ 
tively.  The  low  level  of  NO*  formation  in  flameless  combustion  is 


attributed  to  uniform  low  temperature  in  flameless  mode.  In  the 
other  hand,  the  very  high  temperature  of  the  flame  in  traditional 
combustion  is  responsible  for  high  level  of  thermal  NOx  formation. 
Simulated  results  confirm  that  thermal  NO*  formation  is 
eliminated  in  low  temperature  flameless  combustion,  and  NO* 
formation  occurs  via  N20-intermediate  model. 
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Fig.  13.  N0X  formation  in  the  furnace. 


6.  Conclusion 

Computational  fluid  dynamic  investigation  was  performed  to 
analyze  combustion  characteristics  of  biogas  flameless  regime. 
The  two-step  reaction  scheme  was  simulated  with  EDM  and  the 
standard  k-s  formulation  was  applied  due  to  its  accuracy  and 
robustness  for  wide  range  of  turbulent  flows.  Surface-to-surface 
radiation  heat  transfer  was  modeled  by  DO  due  to  its  great  ability 
in  simulation  of  various  thicknesses.  Due  to  elimination  of  ignition 
in  biogas  flameless  mode,  the  required  enthalpy  for  biogas  auto¬ 
ignition  temperature  is  supplied  by  enthalpy  of  preheated  oxidizer. 
The  effects  of  preheated  temperature  on  conventional  combustion 
were  discussed.  It  was  found  that  in  preheated  biogas  conventional 
combustion,  fuel  consumption  decreases  however  NO*  formation 
increases  drastically.  Indeed,  the  differences  between  reactants 
and  products  temperature  could  intensify  irreversibility  in  tradi¬ 
tional  combustion.  Since  the  temperature  inside  the  chamber  as 
well  as  wall  temperature  is  uniform  in  biogas  flameless  mode  exer- 
gy  loss  decreases  in  flameless  combustion  technique.  Indeed,  NOx 
formation  reduces  in  flameless  mode  due  to  elimination  of  hot 
spots  and  low  level  of  oxygen.  The  streamline  of  the  species  in  both 
conventional  and  flameless  mode  indicate  that  the  species  of  bio¬ 
gas  are  pulled  to  the  oxidizer  streamline  due  to  various  pressures. 
The  mixture  of  the  fuel  and  oxidizer  species  in  flameless  regime  is 
done  faster  than  conventional  mode  due  to  the  turbulence  environ¬ 
ment  and  complete  biogas  combustion  is  projected  in  the  zone 
close  to  the  burner.  Also,  low  concentration  of  02  in  the  flameless 
mode  is  attributed  to  dilution  of  oxidizer  and  the  complete  and 
quick  mixing  process.  Indeed,  very  high  concentration  of  C02  spe¬ 
cies  in  biogas  flameless  products  causes  higher  heat  capacity  and 
better  radiation  heat  transfer  in  the  system. 
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